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Abstract 

The production of c and b quarks in 77 collisions is studied with the L3 de- 
tector at LEP with 410 pb _1 of data, collected at centre-of-mass energies from 189 
GeV to 202 GeV. Hadronic final states containing c and b quarks are identified 
by detecting electrons or muons from their semileptonic decays. The cross sections 
cr(e + e~ — > e + e~ccX) and <r(e + e~ — > e + e~bbX) are measured and compared to next- 
to-leading order perturbative QCD calculations. The cross section of b production 
is measured in 77 collisions for the first time. It is in excess of the QCD prediction 
by a factor of three. 



Submitted to Phys. Lett. B 



1 Introduction 



The measurement of heavy flavour production in two-photon collisions provides a reliable test 
of perturbative QCD because of the large physical scale set by the charm or beauty quark 
mass. Many experiments have studied charm production in 77 collisions |UI21- Beauty quark 
production in 77 collisions is expected to be suppressed by more than two orders of magnitude 
owing to its smaller electric charge and larger mass. At LEP energies, the direct and single 
resolved processes, shown in Figure ^ are predicted to give comparable contributions to the 
heavy flavour production cross section |3J, whereas at low energies the direct process dominates. 
The main contribution to the resolved photon cross section is the photon-gluon fusion process 
7g — > cc(bb). The production rate of c and b quarks in two-photon collisions therefore depends 
on their mass and the gluon density in the photon. Further contributions to charm production 
arise from the Vector Dominance Model (VDM) and from doubly resolved processes and are 
expected to be small. 

This letter describes the measurements of the inclusive open charm and beauty production 
performed with the L3 detector |I] with 410 pb _1 of data, collected at centre-of-mass energies 
from 189 GeV to 202 GeV. 

We identify b and c quarks by tagging electrons or muons from their semileptonic decays. 
Due to the higher mass of the b quarks, these leptons are characterized by a higher momentum 
and higher transverse momentum than those from c quarks. The bb production cross section 
is hence measured using the muon and electron transverse momentum spectra for both leptons 
having a momentum greater than 2 GeV. The charm cross section measurement is made by 
analysing the electron spectrum with momentum down to 0.6 GeV. Muons are always required 
to have a momentum greater than 2 GeV necessary to penetrate the calorimeters and reach 
the muon chambers. 

2 Monte Carlo 

The PYTHIA [3] Monte Carlo is used to model the two photon processes. The non-b events 
are generated with massless matrix elements [H] while for b events massive matrix elements 
are used. The resolved process uses the SaSld photon structure function [Zj. The two-photon 
luminosity function is implemented in the equivalent photon approximation (EPA) [E] with a 
cutoff Q 2 < m 2 . 

Background sources are e + e~ — > e + e~r + r~, e + e~ — > Z/7 — > qq, e + e~ — > t + t~ and e + e~ — > 
W+W~. These processes are generated with JAMVG 0, PYTHIA, KORALZ [TU| and KO- 
RALW ^l] respectively. The detector simulation is performed using the GEANT [12] and 
GHEISHA ^3] packages. The Monte Carlo events are reconstructed in the same way as the 
data. Time dependent detector inefficiencies, as monitored during the data taking period, are 
also simulated. 

3 Event Selection Procedure 

The event selection is performed in two steps. The first one selects hadronic final states pro- 
duced in two-photon collisions, the second identifies a c or b quark by its semileptonic decay. 

^ Electron stands for electron or positron throughout this paper. 
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3.1 Hadronic Two-Photon Events 



Hadronic two-photon events are selected by requiring at least five tracks and a visible energy, 
E vis , below y/s/3. The visible mass, W vis , of the event is calculated from the four-momentum 
vectors of the measured particles, tracks and calorimetric clusters including those from the small 
angle luminosity monitor. These particles are considered to be pions except for unmatched 
electromagnetic clusters considered as photons. W v i S has to be greater than 3 GeV. Figures 
2a and 2b show the E vis and W v i s distributions for the data compared with the Monte Carlo. 
The cut on E vis separates the two-photon process from annihilation processes, characterized 
by high visible energy. The background from e + e~ — > e + e~r + r~ and e + e~ — > t + t~ events is 
suppressed by the requirement on the number of tracks. 

The analysis is limited to untagged events with small photon virtuality. Events are excluded 
when the most energetic cluster in the small angle calorimeters has an energy greater than 0.2 
y/s. Thus the interacting photons are quasi-real: (Q 2 ) = 0.015 GeV 2 , where — Q 2 is the mass 
squared of the virtual photon. 

3.2 Lepton Selection 

Electrons from charm and beauty semileptonic decays are identified by requiring electromag- 
netic clusters in the polar angle range |cos 9\ < 0.725 with momentum greater than 0.6 GeV. 
They should satisfy the following criteria: 

• The candidate cluster is required to match to a track. The difference between the az- 
imuthal angle between the shower barycentre and the track impact point at the calorime- 
ter must be less than 20 mrad. 

• The distribution of energies measured in the crystals of the calorimeter should be com- 
patible with that of an electromagnetic cluster. 

• The E t /pt ratio must be equal to one within 2a, where E t is the projection of the energy 
of the cluster on the plane transverse to the beam, p t is the transverse momentum of the 
track and a is the resolution on this ratio. 

• Photon conversions are suppressed by requiring the distance of closest approach of the 
track to the mean e + e~ collision point in the transverse plane to be less than 0.5 mm and 
the invariant mass of the electron candidate and of the closest track to be greater than 
0.1 GeV. 

After these cuts, 2434 events remain. The background from annihilation processes and 
two-photon production of tau pairs is estimated to be 0.75%. 

Muon candidates are selected from tracks in the muon chambers in the angular range 
|cos 6\ < 0.8. The momentum of the muons must be greater than 2.0 GeV. To suppress 
background from annihilation processes, the muon momentum must be less than 0.1 y/s. After 
all cuts are applied, 269 events remain. The estimated background from annihilation processes 
and two-photon production of tau pairs is 6.0%. 

4 Total cross section cr(e + e _ — ► e + e _ bbX) 

The b cross section is derived from a fit to the data distributions for events where the momentum 
of the muon or electron candidate is greater than 2.0 GeV, corresponding to 269 and 137 events 
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respectively. The b selection efficiency is 2.20% for muons and 1.25% for electrons. The 
transverse momentum of the lepton with respect to the nearest jet, P t , is chosen as the fit 
variable since it has the highest sensitivity to the b fraction. The jets are reconstructed using 
the JADE algorithm ^1] with y cut = 0.1. The energy of the muon or electron is not included 
in the jet. 

A x 2 fit is performed to the data distributions using the sum of the Monte Carlo distributions 
of the non two-photon background Abkg, light quark events iV u ds, c quark events N c5 and b quark 
events N h ^. A three parameter fit is applied in which the number of beauty, charm and light 
quarks are free parameters, whereas A^kg is held fixed according to the Monte Carlo prediction. 
The results of the fit are given in Table ^ which shows also the charm cross section estimate 
from the fit of the muon and electron P t spectra. 

The fit for muons yields a b fraction of 51.6 ± 9.8 (stat) %. As for the electrons, the b 
fraction is 42.3 ± 11.4 (stat) %. The x 2 P er degree of freedom for the muon and electron fits 
are 6.2/6 and 10.1/6 respectively. If no b contribution is included in the fit, confidence levels 
of 2.2 x 10~ 5 and 1.2 x 10~ 3 are obtained for muon and electron respectively. The signal events 
are produced in two separate samples for direct and resolved processes assuming a 1:1 ratio [3]. 
The fitted distributions are shown in Figures 01 and |U 

The resulting cross sections for the luminosity averaged centre-of-mass energy (\/s) = 194 
GeV are: 

a(e + e~ -> e + e-bbX) muons = 14.9 ± 2.8 (stat) ± 2.6 (syst) pb 

a(e + e" -> e + e"bbX) electrons = 10.9 ± 2.9 (stat) ± 2.0 (syst) pb. 

The systematic uncertainties arise from the event selection, jet reconstruction, massive or mass- 
less charm quarks in the event generation, b semileptonic branching ratio, trigger efficiency, 
Monte Carlo statistics and direct to resolved process ratio. Table El shows the values of sys- 
tematic uncertainties due to these contributions. The dominant uncertainty comes from the 
event selection. It is estimated by variation of the cuts and includes detector resolution un- 
certainties and the agreement between data and Monte Carlo. The systematic uncertainty on 
jet reconstruction is assigned by variation of the y cut parameter. The contribution due to the 
uncertainty on the ratio of the direct to the resolved processes is estimated by changing it from 
1:1 to 1:2 or 2:1. 

The combination of muon and electron results gives 

cr(e + e~ — > e + e _ bbX) com ^ nec j = 13.1 ± 2.0 (stat) ± 2.4 (syst) pb. 



5 Total cross section cr(e + e — > e + e ccX) 

In order to increase the statistical accuracy for the charm cross section measurement, the 
electron momentum cut is relaxed from 2 GeV to 0.6 GeV. The cross section is calculated as 
for the data collected at centre-of mass energies from 91 GeV to 183 GeV |2j. In addition, the 
beauty contribution to the number of observed events is subtracted. The open charm cross 
section is then calculated from the number of events with leptons using the equation: 

/Arlept _ yvrlept _ Arlepts 

a(e e — ► e e ccX) = — - — , (1) 

e trig C c 

where the variables are defined as follows: 
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iV tJf is the number of events in the data after the final electron selection (2434 events). 

e tr i g is the trigger efficiency (94.4%) which is determined from the data using a set of 
independent triggers. 

iV b kg is the number of background events which do not originate from two-photon hadronic 
interactions estimated from Monte Carlo (18.3 ± 2.0 events). 

N^ pt is the number of beauty events (169.5) estimated from the cross section measured 
above. 



• C is the total integrated luminosity (410 pb 1 ). 

The c selection efficiency, e' c , is the fraction of c events selected relative to those generated in 
the full phase space. In order to be independent of the Monte Carlo flavour composition, the 
charm purity is written as: 

tt c = (!-—)/(!-— ), (2) 

^data 

where e c and e u d s are the fractions of charm events, Nl ept , and light quark events, accepted 
by the final selection. Equation (2) follows from expressing the number of non-beauty hadronic 
events as: 

Adept — I — Tulept Arlept jylcpt 
c 1 uds c i uds 



+ (3) 

^data ^uds 

where the quantity edata is defined by the following relation: 

Arlept , AT-lept jylept _ Arlept ^lept 

_ iV c + iV uds _ iV obs iV bkg iV b / ,s 

data AThad i AThad Arhad AThad AThad " V / 

iV c ' iV uds JV obs — JV bkg — iV b 

iV b iad , N^ ad and are the number of hadronic events with b-quarks, c-quarks and light quarks 
respectively. and are the observed number of hadronic events and the background 

expectations respectively. This method is insensitive to the absolute normalization of the c 
and uds background Monte Carlo samples, but still depends on the ratio of direct to resolved 
processes in the signal Monte Carlo. 

For the electron sample, the charm purity is 75.0% and the charm selection efficiency is 
0.41%. The charm production cross section in 77 collisions at = 194 GeV is then: 

a(e + e" -> e + e - C cX) electrons = 1072 ± 33 (stat) ± 126 (syst) pb. (5) 

This charm cross section is compatible with the fit results reported in Table H The system- 
atic uncertainties arise from the event selection, direct to resolved process ratio, c semileptonic 
branching ratio, massive or massless charm quarks in the event generation, experimental un- 
certainties on the beauty cross section, trigger efficiency, Monte Carlo statistics and uds back- 
ground estimate. The average charm semileptonic branching ratio used in the simulation is 
0.098 [in]. Table El shows the values of the systematic uncertainties due to different contribu- 
tions. The dominant systematic uncertainty is from event selection and from the variation of 
the ratio of direct to resolved processes. They are estimated as in the beauty study. 

The fit result for beauty production is checked by this counting method in the electron 
case, fixing the charm cross section to the value of equation (5). In addition to the momentum 
cut of 2 GeV, the transverse momentum P t is required to be greater than 1.0 GeV. After all 
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cuts are applied 106 electron candidates remain. The beauty purity is 49.0%, and the selection 
efficiency is 1.2%. The cross section at (y/s) = 194 GeV is 

<r(e + e~ — > e+e ~bbX) e j ec |- rons = 11.3 ± 2.3 (stat) pb, 

in good agreement with the fit result. 

Figure El shows the momentum distribution of the electron candidates. In this plot the 
charm and beauty cross sections predicted by the PYTHIA Monte Carlo are scaled to the 
measured values. The data and Monte Carlo shapes of the momentum distribution show a 
good agreement. 

In the case of muons, where the 2 GeV cut can not be relaxed, the estimation of charm 
production is that derived by a simultaneous fit to the b and c fractions described in the 
previous section: 

a(e + e - e + e"ccX) muons = 814 ± 164 (stat) ± 200 (syst) pb 

at (y/s) = 194 GeV. The efficiency is much lower for the muon sample, about 0.04%, due to 
the higher momentum cut. The dominant systematic uncertainties are the event selection and 
the direct to resolved process ratio. 

The combined value for the open charm cross section at (y/s) = 194 GeV is: 

cx(e + e~ — > e+e_c cX) com ] j nec j = 1016 ± 30 (stat) ± 120 (syst) pb. 

6 Comparisons with QCD Predictions 

The cross sections for open beauty and charm production are compared in Figure 6 to pertur- 
bative next-to- leading order QCD calculations [Sj. The dashed line corresponds to the direct 
process, NLO QCD calculations, while the solid line shows the prediction for the sum of direct 
and resolved processes. The direct process depends upon the QCD coupling constant and the 
heavy-quark mass. The prediction for open charm is calculated using a charm mass of either 
1.3 GeV or 1.7 GeV and the open charm threshold energy is set to 3.8 GeV. The theory pre- 
diction for the resolved process is calculated with the GRV parton density function [TEj. The 
renormalization and factorization scales are chosen to be the heavy quark mass. The direct 
process 77 — > cc is insufficient to describe the data, even if real and virtual gluon corrections 
are included. The data therefore require a significant gluon content in the photon. 

The prediction for open beauty is calculated for a b quark mass of 4.5 GeV or 5.0 GeV and 
the open beauty threshold energy is set to 10.6 GeV. For (y/s) = 194 GeV and a b quark mass 
of 4.5 GeV, this cross section is 4.4 pb. The bb cross section is measured in 77 collisions for 
the first time and is a factor of 3 and about 4 statistical uncertainty standard deviations higher 
than expected. 
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Fit results 


Muon Tag 


Electron Tag 


N bkg 
N bb 
N CE 
N uds 
X 2 / d.o.f. 


16.2 (fixed) 
126.7 ±24.1 
119.0 ±24.0 
n n+ 33 

u - u - 0.0 

6.2 / 6 


2.9 (fixed) 
52.5 ± 14.1 
71.5 ± 14.8 
n n+ 7 - 4 

10.1 / 6 


cr(e + e — > e + e bbX), (pb) 
cr(e + e~ — > e + e~ccX), (pb) 


14.9 ± 2.8 
814 ± 164 


10.9 ± 2.9 
1092 ± 226 



Table 1: Fit to the distribution of the transverse momentum of the lepton with respect to the 
nearest jet. The fit parameters are constrained to be positive. The correlation between N bb - 
and N cE is 75%. 



Source of uncertainty 


Muon Tag 


Electron Tag 


Aa(e + e~ -> e+e'bbX), % 


Aa{e + e~ -> e+e'bbX), % 


Event selection 


14.6 


15.8 


Jet reconstruction 


8.2 


8.2 


Massive/massless charm 


3.0 


3.0 


B(b — > e, /i) 


2.0 


2.0 


Trigger efficiency 


2.0 


2.0 


Monte Carlo statistics 


1.4 


1.8 


Direct / resolved ratio 


1.0 


0.9 


Total 


17.3 


18.4 



Table 2: Systematic uncertainties on <r(e + e — > e + e bbX). 
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Source of uncertainty 


Electron Tag 


Muon Tag 


Aa(e + e' -> e+e-ccX), % 


A ( j(e+e- -> e+e-ccX), % 


Event selection 


8.5 


18.6 


Direct / resolved ratio 


5.7 


10.9 


B(c — > e, ji) 


3.0 


5.0 


Massive/massless charm 


3.0 


3.0 


b background 


2.4 




Trigger efficiency 


2.0 


2.0 


Monte Carlo statistics 


1.9 


6.0 


uds background 


1.1 




Jet reconstruction 




8.2 


Total 


11.7 


24.6 



Table 3: Systematic uncertainties on <r(e + e — > e + e ccX). 
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Direct Single Resolved 

Figure 1: Diagrams contributing to charm and beauty production in 77 collisions at LEP. 
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W„ ifi (GeV) 



Figure 2: a) Total visible energy after applying hadronic selection cuts. A cut of E v i S < 0.33^ 
indicated by the arrow removes most of the backgrounds, b) The visible mass after applying 
hadronic selection cuts. 
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Figure 3: The distribution of the transverse momentum, P t , of the muon candidate with respect 
to the closest jet. 
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Figure 4: The distribution of the transverse momentum, P t , of the electron candidate with 
respect to the closest jet. 
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Electron Momentum (GeV) 



Figure 5: The momentum of the electron candidates. The dotted, dashed and solid histograms 
are the contributions of uds, udsc and udscb quarks from the PYTHIA Monte Carlo. The c 
and b fraction of PYTHIA are scaled to the measured cross sections. 
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Vs (GeV) 

Figure 6: The open charm and beauty production cross section in two-photon collisions. The 
L3 data from both electron and muon events are combined. The statistical and systematic 
uncertainties are added in quadrature. The dashed line corresponds to the direct process 
contribution and the solid line represents the NLO QCD prediction for the sum of the direct 
and resolved processes. 
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